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Abstract 



The effect of a finite nuclear mass on the magnetic moment of the electron 
bound in the ground state of a hydrogen-hke ion is analyzed. Using the 
exact in Za expression for the recoil shift of the energy, I calculate the order 
{Za)^m/M correction to the electron gyromagnetic ratio. The theoretical 
uncertainty in the prediction for the experimentally interesting carbon ion is 
decreased to the level of 1.2 x 10~^. This provides a factor of two improvement 
in the precision of the atomic mass of the electron. 
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Given a particle X with the mass m(X), its relative atomic mass Ar{X) is defined |l|] by 

where ^^C stands for a free neutral atom of carbon 12 in its ground state. At present, the 
most precise value of the electron's relative atomic mass 0, 

A^/=(e) = 0.000 548 579 911 1 (12), (2) 

is obtained by measuring the ratio of the cyclotron frequency of a carbon 12 atomic nucleus 
C(Jc(^^C®^) = 6|e-B|/m(^^C^+) to the cyclotron frequency of the electron uJc{e) = \eB\/me in 
the same magnetic field B. 

It was noted by G. Werth that the high precision measurement of the ratio 
ujl/uJcC''^C^~^), where ujl is the Larmor precession frequency, in the hydrogen-like carbon 
ion i2q5+ i^g^^ offer an alternative source of precise information about the electron mass. 
The idea is to extract this mass from the ratio of the cyclotron frequencies uJc(^'^C^^)/uJc{e) 
by equating the experimental value of the gyromagnetic ratio 

to the corresponding theoretical value 5fth(^^C^+). Such an extraction can be done with 
better precision than in provided that both u!l/uJc{^^C^^) and gth(}'^C^~^) has lower 
uncertainties. 

Using the ratio uJci^^C^+)/uJcie) = 0.000 228 627210 33(50) deduced from the result of 
Ref. 0, the GSI-Mainz collaboration has obtained 0] 

9e.pC^C'^) = 2.001 041 596 4 (8) (6) (44). (4) 

The first two numbers in parentheses show the uncertainties due to statistical and possible 
systematic errors, respectively, and the last number is due to the uncertainty in the ratio of 
the cyclotron frequencies. The theoretical prediction for the above gyromagnetic ratio was 
obtained in Ref. [^: 

c/th(^^C^+) = 2.001 041 590 7 (50) (12) (9). (5) 

Along with the free electron value known with very high precision both experimentally 
0] and theoretically 0, this result includes bound state effects which mainly contribute 
to its uncertainties. Thus, the first theoretical uncertainty is an estimate of the order 
(a/7r)^(Za)^ radiative correction, the second one is due to numerical errors in the first 
radiative correction calculated in to all orders in Za, and the third one is an estimate of 
the order (Za;)^me/m(^^C^^) recoil correction [§,0. 

The prediction (|^) does not take into account the result for the order (a/nY^Za)^ 
correction obtained in Refs. which shifts the central value of (^ by —11.3 x 10~^° 
and decreases the first uncertainty to {a/TT)'^{Za)'^ ^ 2 x 10~^^. In the present work I 
demonstrate that the order {Za)^me/m{^'^C^~^) recoil correction equals 
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(6) 



This correction shifts the theoretical prediction (^) by —1.4 x 10^^^. Remaining uncertainty 
due to the recoil effect can be estimated as me/m(-'^^C^+)(Za)^ ln(l/Za;)/7r ^ 0.7 x 10~^^. 
Including into the results of Refs. and the result (|^), we can deduce from (|[) and 
(@) the improved value of the electron's atomic mass: 

A^/^ie) = 0.000 548 579 909 33 (22) (16) (32), (7) 

where two first numbers in parentheses are the statistical and systematic uncertainties, 
respectively, imported from (1), and the third one is due to numerical errors of the theoretical 
result from Ref. 0. Note that only when added linearly, the uncertainties in (^ and (^ 
cover the difference between the corresponding central values. 

The rest of this paper is devoted to the derivation of the 0{{Za)'^m/M) recoil correction 
(^) arising due to the finite electron-to-nucleus mass ratio m/M. 

To the first order in m/M, dynamics of a hydrogen- like ion is governed by the Hamilto- 
nian 

H = cx{p, - eA{r,)) + f5m - —^+— ' + / dx . 8 

|re — r^vl 2M J 2 

Here a and (3 are the Dirac matrices acting on the spin degrees of freedom of the electron 
with the mass m and the charge e < 0, while Vf. and are the electron position and 
momentum operators, respectively. Corresponding operators for the spinless nucleus with 
the mass M and the charge Z\e\ are and pjq. The electron interacts with the structureless 
nucleus through the instantaneous Coulomb potential. Both particles interact also with the 
transverse vector potential 

A{x) = ?-^ + A{x), (9) 

which consists of the classical and quantum parts. Besides the interaction terms, the latter 
enters into (|) through S{x) = —i5/5A{x) and B{x) = —i'Vx x A.{x). 

To find the order m/M correction to the electron magnetic moment, we have to split 
the Hamiltonian into an unperturbed part and a perturbation. The former describes the 
system in the limit of no recoil, M — > cxd, and hence depends on the position of the electron 
with respect to the nucleus r = r^. — r^. Since in this limit the nuclear momentum 
drops out of the Hamiltonian (^), the canonical momentum of the electron p^ coincides 
with p = —id/ dr. The most appropriate way to maintain the dependence on r in the 
unperturbed Hamiltonian is to perform the unitary transformation [|1^] 

ij^Uij, U = exp (^ieB ^^^^^ ^ . (10) 

This transformation simplifies the perturbation theory which otherwise should include the 
virtual transitions into excited P states of both the electron with respect to the nucleus and 
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the ion as a whole with respect to the origin. In fact, to the first order in the magnetic field 
(sufficient for the analysis of the magnetic moment), 

U^l+ leB'^^i^ = 1 + leB?^, (11) 

where the vector R points to the center of mass. Thus we see that the unitary transformation 
(p!0|) describes the admixture of P states induced by the external magnetic field, in both the 
relative and the center of mass motions. 

The transformed Hamiltonian naturally splits into the unperturbed part and the pertur- 
bations: 

U^HU = Hb + Kad + Kec; (12) 

/ Bxr\ ^ Za f , £^ + B"^ 
HB = a\p-e — - — \ + pm \- J dx , (13) 

V.ad = -eaA(r), (14) 
_ (p^ + ZeA{0)f {eB xr,p^ + ZeAjO)} 

2M 4M ■ ^ ^ 

The Hamiltonian Hb describes the electron in the static external field (Coulomb plus mag- 
netic) and the free electromagnetic field. In the limit B ^ it turns into the Dirac Hamil- 
tonian in the Coulomb field, Hq. The perturbation V^ad induces the electron interaction 
with the quantized electromagnetic field, while the perturbation Vrcc gives rise to the recoil 
correction. We can put rjy = in Vj-ec which already contains the overall 1/M. Looking 
for the linear response of the system to the external magnetic field, we can also neglect the 
order B^ term in V^ec- 

At zero magnetic field, the recoil correction to the energy of the Dirac electron in the 
Coulomb field was expressed in terms of the solution to the Dirac-Coulomb problem in Refs. 



2M 



2\ 



B=0 



l^iHiP- D^) Geo+. (p - D^) l^o) . (16) 
J Ztxi 



1 I' doj 



The average value on the left hand side is calculated over the fluctuations of the quantized 
electromagnetic field, with V^^d included an appropriate number of times (see [|ll|-|T^ for 
details). The average value on the right hand side is calculated over the eigenstate of the 
Dirac equation in the Coulomb field at -B = 0, 

(Z cx \ 
cxp + pm--^JijQ{r) = EaMr), (17) 

the operator Z)^ describes the transverse (magnetic) exchange between the electron and the 
nucleus, 

^ / X ^ f dk Anexpiikr) ( k(a.k)\ 
and is the Green function for the Dirac equation in the Coulomb field. 
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(^E-ap-(3m+ ^) ^^(r, r') = 6{r - r'). (19) 

The contour of integration over uj in (|16D goes from minus infinity to zero below the real 
axis, rounds zero from above, and then proceeds to plus infinity above the real axis. 

In Ref. [0 , Eq. ([1^) is derived from the requirement of the recoil correction's invariance 
with respect to gauge transformations of the quantized electromagnetic field. That derivation 
does not exploit a particular form of the potential entering into the Dirac equation and hence 
can be directly generalized to the present case of the Coulomb plus magnetic field: 

(Kec) = -j^J^^ (^b\ (7r+ - D^) G%^^^ (77+ - 1^^) . (20) 

The Dirac equations for the wave and Green functions iI)b and differ from (|1^ and 

(p!9|) by the substitution 



In Eq. (pOD, the operator 



eB X r , , 

P — > TT = p . (21) 



tB X r , ^ 

'K^ = p+^— (22) 



arises instead of tt due to the second term in the perturbation operator V^ec, (|T5|). The same 
result as ( |20| ) was obtained recently in Ref. |]T5[ in a different manner. 



An equivalent form of Eq. (|20|), which is more convenient for the further analysis, stems 
from the Dirac equation: 

(V^ec) =(Kec)low+(Kcc)high, (23) 

= If I ^A{-\^^^Hs\ - D^ G|^+^ (-[7T^,Hs] + D^) . (25) 



^^'s^ M J 2m WW - J ~" -- \uj- 

Here and below {■ ■ ■) ^ denotes the average value over ipB- In what follows I demonstrate 



that in perfect analogy with the case of zero magnetic field, (^J) includes the lower order 
{{ZaY and {ZaY) contributions to the energyQ, while the expansion of (|25D starts with 
{Zaf. 

The average value of the local operator ( p^ can be easily found in an analytic form. 
Taking the square of the operators in both sides of the Dirac equation, 

Zoi \ 

a(7r+ -eBxr)- -yj ijjB{r) = {Eb - (3m) ^/'^(r), (26) 



^It concerns the S states; for L > 0, there is also a purely kinematic 0{1) effect |16]. 
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and taking also into account that 



Za ' 



r J 



+ Za 



a. X n 



2r2 



= r X 7r_i 



(27) 



we get 



low 



2M 

1 / ex. X TL 

+ — ( eBj + ZaeB{oL x n) - Za—^L 



M 



2r2 



[Za) 
2r2 



(28) 



B 



where j = 1 + 12/2, while I = r x p. The sum of the last two terms in (|28|) vanishes, 

(Za?\ 



-Za^L 
2r2 



2r2 4 



{[HB,Za{n+,n}])^ = 0, 



so that 



(^rec) 



El + m^ - 2EBm {p)^ + 2eBj + 2ZaeB {a x n) 



B 



low 



2M 



To the first order in the magnetic field, 



Eb Eq — go 
{P)b - 



eBj 



2m ' 
dHs \ OEb 



dm 



B 



dm m 



eBf 

2m 



ZaeB {a x n) ^ ZaeB {a x n)^ 



AeBj (El 



(29) 

(30) 

(31) 
(32) 

(33) 



Here gn is the bound electron gyromagnetic ratio without radiative and recoil corrections. 
The average value in (^3]) is calculated for the electron S states by using virial relations 



for the Dirac equation ||T^. In the ground state, go = 2/3(1 + 2wl — (Za)^) [jT8| and 



Eq = my 1 — (Za)2 (see, e.g., [|T9l) so that the recoil correction to the gyromagnetic ratio 
derived from (123) reads 



2m 
3M 



1 - (Za) 



l + 2{Za) 



(34) 



and agrees with the result obtained in Ref. []15|. The first term of its expansion in Za, 



5i^)g = !^(Zaf, 



(35) 



reproduces the result of Refs. p0|JT0|j21 



In order to prove that the expansion of (|25|) starts with {ZaY, let us first note that 

Zan 



(36) 
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and hence the only difference between (^5]) and its S = counterpart is due to first order 
differences between ipB and ipQ and between and G. For the hard scale contribution to 
(p5|), when a; ~ m, the interaction of the highly excited electron with the magnetic field 
described by G^ — G, is an order (Za)^ effect. In fact, since at the hard scale Bxr \B\/m, 
the inclusion of this interaction does not change the power counting valid for the case of the 
zero magnetic field. 

One could suppose that the B = power counting breaks down at the atomic scale 
where r ~ l/[mZa). The following arguments, however, show that this is not the case. If 
Za <^ 1 we can treat atomic scale effects in the nonrelativistic approach, i.e. expanding 
in p/m ~ Za. In particular, such an expansion for the effective potential describing the 
nonrelativistic electron interacting with the magnetic field, starts with the Pauli term, 

where g^ree = 2 + a/ir + . . . is the gyromagnetic ratio for the free electron and cr/2 is the 
nonrelativistic spin operator. The perturbation ( P7D is diagonal in the basis of solutions to 
the Schrodinger equation in the Coulomb field. Therefore, it has no effect on those solutions 
as well as on the corresponding Green function and hence on the atomic scale contribution 
to ®. 

The next term in the expansion of the effective potential, of order /m? ~ [ZaY, was 
found in [||,|: 

Here E = —Zen/r"^ is the Coulomb electric field and = p^/(2m) — Za/r is the 
Schrodinger Hamiltonian in the Coulomb field. As far as V^"^^ ~ {Zay\B\/m and the charac- 
teristic energy denominator is \Eo—En\ ~ m(Za)^, this perturbation induces 0{\B\/m'^) cor- 
rections to the Schrodinger wave function and to the corresponding Green function. Hence, 
the first term of the expansion in Za has the same order (Za)^ for both and |-B|^ 
contributions to (|25|). 

Thus we come to the conclusion that the order (Za)'^ recoil correction to the gyromag- 
netic ratio of the bound electron can be drawn from the expansion of which results in 
Eq. i). 

The next step in the reduction of uncertainty of the electron mass (|^) can be made by 
a more accurate evaluation of the first radiative correction to the g factor of the bound 
electron. As far as Z = 6 is yet much less than 1/a 137, an analytic calculation of the 
order (a/7r)(Za)^ correction can also prove to be quite useful. 

In conclusion, the new value of the electron's atomic mass is extracted from the recently 
measured magnetic moment of the electron bound in the hydrogen-like carbon ion. Due 
to the high precision of this measurement, as well as of the theoretical prediction for the 
gyromagnetic ratio of the bound electron, the overall uncertainty of the new value is reduced 
by a factor of about two as compared with the previous result extracted from the ratio of 
the carbon nucleus cyclotron frequency to that of the electron. 
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